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Abstract 

We present a systematic approach to the Unearised Yang-Mills-Higgs equations in 
the background of a 't Hooft-Polyakov monopole and use it to unify and extend 
previous studies of their spectral properties. We show that a quaternionic formula- 
tion allows for a compact and efficient treatment of the linearised equations in the 
BPS limit of vanishing Higgs self-coupling, and use it to study both scattering and 
bound states. We focus on the sector of vanishing generalised angular momentum 
and analyse it numerically, putting zero-energy bound states. Coulomb bound states 
and infinitely many Feshbach resonances into a coherent picture. We also consider 
the linearised Yang-Mills-Higgs equations with non-vanishing Higgs self-coupling 
and confirm the occurrence of Feshbach resonances in this situation. 

1 Introduction 

Solitons - spatially localised and stable solutions of non-linear differential equations - are 
widely used in modelling physical and biological phenomena. In many of these appli- 
cations, it is also of interest to study the behaviour of small perturbations around the 
soliton. In the context of quantum field theory, the small perturbations are interpreted 
as particles after quantisation [1]. In that case, studying perturbations around a soli- 
ton amounts to studying the interaction of these (quantum) particles with the classical 
soliton. In Skyrme's model for nuclei [21 El S], for example, the solitons describe nuclei 
and the small perturbations correspond, after quantisation, to pions. Their properties 
in the background of a Skyrmion can therefore be interpreted in terms of pion-nucleus 
scattering. 

Mathematically, studying perturbations of a soliton amounts to studying the spectral 
properties of the linear operator obtained by linearising the soliton equations around the 
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soliton. It turns out that these spectral properties are often interesting, see [T] for a 
textbook treatment of some examples. A particularly interesting example is provided by 
magnetic monopoles in SU{2) Yang-Mills-Higgs (YMH) theory. It has long been known 
that the linearised YMH equations in the background of a monopole have interesting zero- 
modes in the BPS limit [HI El [7] and that they support infinitely many bound states [H]. 
Scattering of fermions off the monopole has also been studied extensively P HO], HU [12] . 
However, fairly recently interesting new scattering phenomena were observed. In [13] and, 
more recently and in more detail in [13], Fodor and Racz studied spherically symmetric 
but non-linear perturbations of the 't Hooft-Polyakov monopole, in the BPS limit. The 
authors studied the evolution of such excitations numerically and found that the monopole 
holds on to a significant fraction of the energy from the initial excitation for a surprisingly 
long time. In particular the amplitude of the excitation decays as t~^^^ for late times. 

Motivated by this work, Forgacs and Volkov carried out a perturbative analysis of the 
't Hooft-Polyakov monopole, still in the BPS limit and preserving the hedgehog form of 
the monopole. Their linearised equations of motion are a system of two weakly coupled, 
second order ordinary differential equations. By using the dispersion relation, one of 
these channels is seen to be massive and one massless. Some intuition about the system 
is gained by artificially decoupling it and considering only the massive channel. The 
potential appearing in it has an (attractive) Coulomb tail, so that it possesses infinitely 
many bound states approaching a critical value, beyond which there is a continuous 
spectrum. On re-coupling and considering the full system, a phase shift analysis of the 
massless channel shows that the infinity of bound states in the massive channel turn into 
an infinity of resonances in the coupled system. The energy held by the resonances leaks 
out slowly to the massless channel, leading to the long-lived excitation observed by Fodor 
and Racz. 

The aim of this work is to extend the results of Forgacs and Volkov and to place it in 
the context of other spectral properties of the hnearised YMH equations. We develop a 
quaternionic formalism for studying the linearised YMH equation in the case of vanishing 
Higgs self-coupling (i.e. in the BPS limit). We show that this quaternionic formalism 
allows for a systematic treatment of perturbations, organised in terms of the eigenvalues 
of the generalised angular momentum operator (combining orbital angular momentum, 
spin and isospin). For vanishing generalised angular momentum we recover the equations 
studied by Forgacs and Volkov but also find another system, consisting of two coupled, 
second order ordinary differential equations. These have bound states, already found by 
Bais and Troost in [8]. The two systems found for vanishing generalised angular mo- 
mentum thus already display a wealth of interesting spectral properties, including bound 
states at zero energy, bound states embedded in the continuum and resonance scattering. 
Moreover, the BPS condition allows one to map either of the systems into equivalent but 
sometimes simpler systems using a first order differential operator (essentially a super- 
symmetry charge, but we do not consider the fully supersymmetric theory here). This 
turns out to explain some of the surprising features we find. 
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We also consider the case of non- vanishing Higgs self-coup hng A. In this regime, our 
quaternionic formalism is no longer effective. Thus we do not study general perturbations 
but instead focus on the generalisation of the system studied by Forgacs and Volkov 
when A 7^ 0. We find that the Coulomb tail in the massive channel is replaced by an 
attractive potential. This is strong enough to support infinitely many bound states 
after decoupling (by hand), and our numerical analysis suggests that it will also produce 
infinitely many resonances in the coupled systems. We conclude that the qualitative 
features found by Forgacs and Volkov survive the "switching on" of A. 

The paper is organised as follows. In Sect. 2, we study the general form of the linearised 
YMH equation around a background which satisfies the first order BPS equation (which 
includes the 't Hooft-Polyakov monopole in the BPS limit). We introduce a quaternionic 
language and show that the linearised YMH equations for stationary time- dependent per- 
turbations can be expressed as a quaternionic wave equation, supplemented by a back- 
ground gauge condition. For non-zero kinetic energy there is a second, equivalent form 
of this wave equation, related to the original one by the application of a Dirac-type op- 
erator. Our strategy for studying the linearised YMH equations is therefore to study 
the quaternionic wave equation and to check if solutions satisfy the background gauge 
condition. 

In Sect. 3, we carry out a partial wave analysis of the quaternionic wave equation and 
derive two systems of two second order ordinary differential equations which arise in 
the sector with vanishing generalised angular momentum. We also derive the form of 
the alternative but equivalent systems obtained by acting with the Dirac-type operator 
of the previous section. Sect. 4 contains a detailed, numerical investigation of the two 
systems found in Sect. 3. Even though these systems look superficially very similar, their 
spectral properties are quite different. One of the systems is the hedgehog system already 
discussed in [15]; we briefly repeat the analysis of this system, using it to set out our 
conventions. We then observe that the other system decouples, after application of the 
Dirac-type operator, into one channel which supports bound states and another which 
only has scattering states. We are able to relate the bound states to those discussed in 
[H]. The scattering states do not satisfy the background gauge condition and therefore 
are not valid bosonic states in the linearised theory, but we point out their relation to 
the fermionic scattering states studied in [12]. In Sect. 5, we continue our study of 
SU{2) monopoles, but allow for non-zero Higgs self-coupling. We perturb around the 
background of the non-BPS 't Hooft-Polyakov monopole and find scattering resonances 
in the linearised hedgehog fields. Finally, in Sect. 6 we discuss possible extensions of our 
work and the interpretation of our results in the context of electric-magnetic duality. 



3 



2 Perturbing the BPS Monopole 



2.1 The BPS Monopole 

Much background material for this section can be found in the textbook ||4j, to which we 
refer for details. We work on four-dimensional Minkowski space-time with coordinates 
x^, fi = 0, 1,2,3 and Minkowski metric rj^^, = diag(l, —1, —1, —1). We denote the time 
coordinate by x° or t and write three-component Euclidean vectors with bold letters, 
e.g. X = {xi,X2,X3). We write the inner product as x ■ y and r for the spatial radial 
coordinate, i.e. r = |x| = ^xl + xl + xl The YMH model we are interested in has 
gauge group SU{2), often referred to as isospin symmetry. Monopoles emerge when the 
Higgs mechanism breaks the SU{2) symmetry to U{1). In the notation and nomenclature 
of this paper we treat this U{1) as the gauge group of Maxwell electrodynamics. We work 
in units where the speed of light and the gauge coupling are 1. 

The fields of the YMH model are an SU{2) gauge potential A^, coupled to a Higgs field 
(J). Both take values in the Lie algebra 3u{i) and transform in the adjoint representation 
of SU{2). For the Lie algebra su(2), we use the basis ta = —-^Ta, where Ta,a = 1,2,3, 
are the Pauli matrices, with brackets [ta,tb] = ^abctc, noting that this is not the same 
convention as [4J. We will also need an inner product (, ) on 5u{i), which we normalise so 
that {ta,tb) = 6ab- The covariant derivative is = + [A^, ] and the Yang-Mills field 
strength tensor, or curvature 2-form, is F^^, = d^A^ — d^A^ + [A^, Ay\. From this we can 
extract the non-abelian electric field Ei = Fqi and the non-abelian magnetic field 

Bi = -]^eijkFjk, i,j,k = l,2,3. (2.1) 
The YMH Lagrangian density is 

^ = -\{F,uFn + l{D,<i>Dy) - ^(1 - (2.2) 

where |0p = (0,0). The equations of motion derived from the Lagrangian density are 

D,,Dy = A(l - 101^)0, (2.3a) 
D^F'^'^ = [D>,0]. (2.3b) 



Static configurations play an important role in this paper as background configurations. 
For such configurations we work in the temporal gauge Aq = 0, and assume time- 
independence of the remaining fields Ai and 0. Sometimes we collect the static gauge 
field into a spatial one-form A = Aidx\ and write (A, 0) for the field configuration. For 
such configurations, the energy computed from the Lagrangian (2.2) can only be finite if 
we impose the boundary condition 

101 1 as r oo. (2.4) 
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The equations of motion (2.3) then reduce to 



AA0 = -A(l - 101^)0, 
-[^0,0]. 



DiFij 



(2.5a) 
(2.5b) 



The 't Hooft-Polyakov static monopole solution found in fW\ [T7] has the hedgehog form 



A;(x) 



Xk, 



1 - W{r))ea^kta, 



_ H{r) 

(X) — —p^Xata- 



(2.6) 



Regularity at the origin requires that W{0) = 1 and H{0) = 0, while the boundary 



conditions at infinity (2.4) are satisfied if we require H — r as r — )■ oo. To recover a 



more standard definition of the Higgs field H{r) we can scale by — r. Our definition of 
H{r) follows the conventions of Forgacs and Volkov in [IS] for the hedgehog ansatz, for 
better comparison with their work. 



H(r) and W{r) satisfy differential equations derived from (2.5), 



-r^f^ + W' + H'-l]W = 0, 



(2.7a) 
(2.7b) 



The BPS limit amounts to setting A = in (2.2) but maintaining the boundary condition 



on 0. As the mass of the Higgs is proportional to a/A, the Higgs field is massless in this 



limit. The static field equations (2.5) become 



AA0 = 0, 

DiFij ■ 



while the equations (2.7) reduce to 



-r^^^ + VT^ + if^-l 1 = 0, 



— r 



dr"^ 



+ 2W^ \H = Q. 



These equations have an analytic solution found by Prasad and Sommerfield in [18 

H{r) = 1 - r coth(r), W{r) 



(2.8a) 
(2.8b) 



(2.9a) 
(2.9b) 

% 

(2.10) 



sinh(r) 

The corresponding field configuration is called the Bogomol'nyi-Prasad-Sommerfield (BPS) 
monopole. 

Further insight into the BPS limit can be gained from considering the static energy 

1 



E=-l {Bi, Bi) + (A0, A0)rf'a;, 



(2.11) 
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and rearranging it to 

E=^J{Bi + A0, Bi + D4)d'x - j di{Bi, (2.12) 

where we used the Bianchi identity DiBi = 0. Using Stokes' law and the quantisation of 
magnetic flux, one flnds 



E=^j{B, + D4, Bi + D4)d^x + 27rAr, 



(2.13) 



where is an integer called the monopole number. Then for > (i.e. monopoles as 
opposed to anti- monopoles), the energy is bounded by 



E > 2tiN. 

This bound is saturated when the Bogomol'nyi equation [19 

B, + = 



(2.14) 



(2.15) 



holds. The exact monopole solution (2.10) satisfles this equation, and W{r) and H{r) 
thus satisfy flrst order equations derived from ( |2.15 ), which we will use later: 



rH'ir) = W{rY + H{r) - 1, 
rW'{r) = W{r)H{r). 



(2.16a) 
(2.16b) 



It is straightforward to check that the first order partial differential equations ( 2.15[ ) imply 



the second order equations (2.8), and that the first order ordinary differential equations 



(2.16) imply the second order equations (2.9). 



2.2 Static linearisation 

We fix a static background configuration (71^,0*), assumed to satisfy the Bogomol'nyi 



equation (2.15), write Df for the covariant derivative Df = di + [Af, ] and likewise 
Ffj for the curvature of Af. Although we are ultimately interested in time-dependent 
perturbations, we begin by considering static perturbations (a, ip) of {A^, 0"*) and insert 



Ai — A^ + Oj, 



(2.17) 



into ( |2.8[ ). Substituting into ( |2.8a[ ), applying the static background gauge condition and 

] = 0. (2.18) 



collecting linear terms in the perturbation we find 



To linearise (2.8b), we need to use that, to linear order, 

Fij ~ F^j + diaj - djai + [ai, A^] + [A^ aj]. 



(2.19) 
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Then substituting into (2.8b), the left hand side is 



Likewise the right hand side of (|2.8b|) becomes 



(2.20) 



(2.21) 



Putting (2.20) and (2.21) together and applying the static equation (2.8b) we obtain the 



linearisation 



(2.22) 



In most of the remainder of the paper we will only use the covariant derivative Df and 
the field strength F^^ associated with a fixed static background configuration. In order to 
simplify notation we therefore drop the superscripts on A^,(f)^, and also on the associated 
covariant derivative and curvature. In this notation, the linearised Yang Mills Higgs 
equations for static fields are 



DiDiip + [flj, A0] + Di[ai, 0] = 0, 
DiDiUj - DiDjai + [a^, Fij] -- 



D,^] + Dj(f)] + [(f), [a 



'J' 



(2.23a) 
(2.23b) 



We can similarly linearise the Bogomol'nyi equations. Substituting (2.17) into the Bo 



gomol'nyi equation (2.15), linearising, using the Bogomol'nyi equation and re-naming 
(j)^ — )■ 0, y4| — )■ Ai gives 

eijkDjak = Diip + [at, (f)]. (2.24) 

There are infinitely many solutions of this equation with = —Did, = [d, </>], where i) 
is an arbitrary function on with values in su{z). These do not physically change the 
original static solution, as they are infinitesimal gauge transformations. We can exclude 
such solutions by requiring that perturbations {ai,(p) satisfy 



(2.25) 



for all 'i? which are non zero on a closed and bounded subset of R^. The requirement 
of compact support means that we can integrate by parts and rearrange to obtain the 
background gauge condition 

A«i + [0, ^] = 0. (2.26) 

Interestingly, the linearised Bogomol'nyi equations together with the background gauge 
condition imply the linearised YMH equations (just as solutions of the Bogomol'nyi equa- 



tion (2.15) are solutions to the static YMH equations (2.8)). To see this, we apply A to 



(2.24) and then use (2.15) to obtain (2.23a). In order to derive (2.23b), we apply Dj to the 



background gauge condition (2.26) and [0, ■] as well as Di to the linearised Bogomol'nyi 



equation (2.24). The algebra is a little tedious, and makes repeated use of (2.15). We 



will give a much quicker derivation in the next section. 
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2.3 Quaternionic Formulation 



We will now show that the language of quaternions is very convenient for studying the 
linearised equations of the previous section. We denote the set of all quaternions as H 
and introduce the usual basis Cq,, a = 1,2,3,4. The real unit quaternion 64 commutes 
with all quaternions and is often written as a 1 or omitted. The remaining (imaginary) 
quaternions satisfy 

eiCj = -6ij + eijkCk {i, j,k = 1,2,3). (2.27) 

We can identify ej = —icj, where aj are again the Pauli matrices (but not denoted Tj here 
in order to avoid confusion with the isospin Lie algebra), and 64 with the 2x2 identity 
matrix I2. The conjugates are 

= — Cj, 64 = 64. (2.28) 



We combine the gauge fields Ai and into a quaternion-valued field 

(l=Aiei + (j). (2.29) 

Since Ai is an isovector-vector and (p is an isovector-scalar, we can view this field as a 
map 

Q,: ^ H (8)SU2. 

Next, we define Dirac-type derivative operators 

Ip = Ae,, + [0, ], IP^ = Ae. - [0, ], (2.30) 

which act on functions ^ : — H ®SU2 by quaternionic multiplication on the quaternions 
and by commutator on the isospin part SU2. These operators are closely related to the 
Bogomol'nyi equation for the background field Q,. Note that 

Jp^lp = -D^ - 02 + ^D,(/) - A)e„ Iplp^ = - A' - 0' + (-A0 - A)e„ (2.31) 

so that, for BPS monopoles, 

Ip^lp = -A - 0' + 2Aei, $1/)^ = -A - 0', (2.32) 



by virtue of (2.15). 

It turns out that the linearised YMH and BPS equations can both be expressed very 
compactly in terms of these operators. To see this, let 

q = aiCi + if. (2.33) 

We look for the quaternionic expression equivalent to the linearised Bogomol'nyi equation 
(2.24) and observe that 

Ipq = Ditti + [0, ip] + ei{Diip - eijkDjak - [0, ai]). (2.34) 
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The real part of (2.34) vanishing is precisely the background gauge condition (2.26) 



while setting the complex part to zero is equivalent to (2.24). The linearised BogomoPnyi 



equation (2.24) togther with the background gauge condition therefore have the following, 



very simple quaternionic formulation 



This was used extensively in 



U)q = Q. 



(2.35) 



We now express the linearised YMH equations (2.23) in quaternionic notation, which has 



not been previously considered. We expect it to be some second order equation in Ip, 
since the linearised Bogomol'nyi equation is a first order equation in ^. In fact, we will 
now show that the quaternionic equation 



WIp'q = 



(2.36) 



is equivalent to the linearised field equations ( |2.23a ) and (2.23b), provided the background 



gauge condition (2.26) holds 



In order to prove our claim we first note a number of useful relations. Using the definition 
of the curvature as the commutator of covariant derivatives, as well as the Bogomol'nyi 
equation (2.15) and the definition ( |2.1 ) of the non-abelian magnetic field, one finds 



eiimDiDiip = ^eiim[Di, Di\Lp = ^eiim[Fii,ip] = [Bm,'f] = [^,Dm(l)]. (2.37) 



The Leibniz rule, again with the Bogomol'nyi equation (2.15) and the definition (2.1) of 
the non-abelian magnetic field, gives 



Then we compute 

Ip^lpq = lp\Diai + [0, + lp^ei{D,ip - eijkDjak - [0, a^]) 
= lp\Diai + [<P,^]) 

- DiDiif - [ai, A0] - Di[ai, 4>] 

+ DiDiajCj - DiDjaiCj + ej[ai, Fij] - ej[(j), Djip] - ej[(p, Djc, 



(2.38) 



-ej[(p, [a 



7' 



(2.39) 



The terms after the last equality sign are, in the first line, the operator Ip = eiDi 



applied to the the background gauge expression (2.26), in the second line the linearised 



field equations (2.23a), and in the third line the linearised field equations (2.23b). Thus 



with the background gauge condition Ditti + [0, (f] = imposed, the quaternionic equation 



(2.36) and the linearised static YMH equations (2.23) are equivalent, as claimed 



In the quaternionic notation it is obvious that the first order equation Ipq = (combin- 
ing background gauge condition and linearised Bogomol'nyi equation) implies the second 
order equation Ip'^ Ipq = (and hence the linearised static field equations, since the back- 
ground gauge is in place). 
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2.4 Time dependent perturbations 

We now introduce time-dependent perturbations around a static configuration {A^, 0^*) 
satisfying (2.15), using the following stationary ansatz 

A{t,^) = 4^(x) + a,(x)e*'^*, 0(t,x) = 0^(x) + y.(x)e*'^*. (2.40) 

We recall that we work in the temporal gauge Aq = and the BPS limit A = 0. Inserting 
the ansatz in the YMH equations (2.3), linearising and re-naming again 0^ — )■ 0, — )■ Ai, 
we find 



DiDiif + [ai, Di(f)] + Di[ai, 0] 
DiDiaj - DiDjUi + [ai, Fij] 



Djip] + [v9, Djcj)] + [0, [aj, 0]] - u'^aj 



(2.41a) 
(2.41b) 



For = we recover the static linearised YMH equation (2.36), as one would expect. 

With the results of the previous two subsections, we can express the linearised equations 
in quaternionic language. Provided the background gauge condition (2.26) holds, the 
following simple equation 

(2.42) 



is equivalent to the stationary linearised YMH equations (2.41). This observation is 
fundamental for the remainder of this paper, and the foundation of our strategy for 
investigating (2.41 ) by studying (2.42 ) and then imposing the background gauge condition. 



Like all differential equations, the equation (2.42) can be written as a first order system. 
In this case, it takes the form of a Dirac equation: 



U) -dt 



(2.43) 



where q and p are quaternion valued fields. 

With the stationary ansatz we have the eigenvalue equations 



U)^ 

Ij) 



(2.44) 



This equation is closely related to a Dirac equation studied by Jackiw and Rebbi in [5], 
where they investigated the zero-energy solutions. Bais and Troost |8] looked at the same 
Dirac equation, extending the analysis to higher energies, and looking for bound states. 



The other second order equation which it is possible to arrange (2.44) into is 



(2.45) 



This equation is equivalent to (2.42), provided tu ^ 0. The map between solutions is 

l/)q = up, and l/)^p = uq, (2.46) 



10 



as seen in (2.44). Explicitly, we can obtain a solution of (2.45) from a solution of (2.42) 



by applying ij) to each side of (2.42) and using (2.46): 



(2.47) 



where we also divided by oj. Conversely, we can also map a solution of ( 2.45[ ) into a 



solution of (2.42) by substituting for p using (2.46), provided tu ^ 0. As we shall see later 



in this paper, it is fruitful to investigate both (2.42) and (2.45), using the map (2.46) to 
relate the results. 



Finally, we note that the transformation (2.46) also provides a convenient way of checking 



if a solution q of (2.42) satisfies the background gauge condition. According to (2.34), 
the latter is the requirement that the real part (in the quaternionic sense) of Ipq vanishes. 
In other words, to see if the background gauge condition holds we simply check if the 



quaternion p obtained from q according to (2.46) has a vanishing real part. 



3 The quaternionic wave equation in the background of the BPS 
monopole 

3.1 Structure and symmetries of the wave equation 



We begin our detailed study of the equation (2.42) in the case where the background field 
is the BPS monopole (2.6), with the profile functions given in (2.10). The Dirac operator 

W), ... H 



introduced in (2.30) now takes the form 

0\ 



"^BPS 



ieidi) + 



e X X 



xt. 



Using the equations ( |2.16[ ) satisfied by the profile functions (2.10) we find 

H 



Xk. 



W)eaikta 



H 



XiXa 



(H 



H\ XkXb 



^aik^abctc 



-r-rl XiXa,^ ^-^r^H WH 

„ H'ta- ^{l + W)-ta+ 



SO that 
{Ip^lp) 



BPS 



-D 



-A- 



2(1 - W) 



L t- 



[l-Wf 



t^ + 



x-t 



ly ^ ry ^ ry ^ 

+ — ■ t)(x ■ e) - ^ M x ■ t)(x ■ e) + — (e ■ t), 
where the orbital angular momentum operator L has components 

Li (^ijkXjOk, 



(3.1) 



(3.2) 



(3.3) 



(3.4) 
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and the Laplace operator can be written as 



A = dl + dl + dl 



We also note that the differential operator in (2.45) now talces the form 



I BPS 



-Di - 0^ 



[l-Wf 



wf 



fx-tl 



(3.5) 



(3.6) 



We will use the quaternionic formulation for studying the operators in (3.3) and (3.6). 
This means that we will let them act on functions 



(3.7) 



as explained after (2.30). We begin with (3.8), which arises directly from the linearisation 



discussed in the previous section, and look for eigenfunctions, i.e. solutions of 

{U)^U))BPsq = UJ^. 



(3. 



The BPS monopole is spherically symmetric in the sense that a spatial rotation can be 
compensated for by an iso-rotation. The operator generating the combined spatial and iso- 
rotations can be expressed in terms of the angular momentum operator L ( 3.4[ ), the spin 
operator s = (whose components act on the quaternion part of (3.7) via commutator), 
and the isospin operator t (whose components ti = —^Ta act in the adjoint representation 
on the 5U2 part of q). It is easy to check that, with our conventions, the components of 
the generalised angular momentum operator 



L + s + t 



(3.9) 



satisfies [Jj, Jj] = eijkJk and commutes with it) bps and Ipsps defined in (3.1). As a result 



we are able to organise eigenfunctions of the operator [ij)^ IJ))bps in terms of multiplets 
of the generalised angular momentum operator J. Doing this in practice is the subject of 
the next section. 



3.2 Partial wave analysis 



In order to split the set of eigenfunctions in (3.8) into irreducible representations (multi- 



plets) of the generalised angular momentum operator J we apply basic results from the 
representation theory of SU{2). Denoting irreducible representations of SU{2) by their 
spin j G |N we recall the basic tensor product decomposition rule 



31® 32= n. (3.10) 

n=lii-i2| 
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From the point of view of representation theory, we may think of quaternions as a direct 
sum of a spin and a spin 1 representation of s: 



Hi ~0©1. 



The quaternionic function (3.7) can be viewed as a tensor product of (i) a scalar function 
on R^, (ii) a quaternion and, {in) a spin 1 representation of the isospin operator t. The 
three terms in J act on each of these separately according to J = L(8)l(8>l+l(8)s®l+l®l®t. 
If we split the space of scalar functions on further into a tensor product of functions 
of the radial coordinate r and the space LP'{S'^) of square-integrable functions on the 
two-sphere, then J acts trivially on the radial functions, and the decomposition of L^(S'^) 
into irreducible representations of the orbital angular momentum operator L is the usual 
decomposition of functions on the two-sphere into spherical harmonics: 

oo 



Thus, applying the rule (3.10) to the tensor product of quaternions and 5U2 we observe 

(0 © 1) (g) 1 = 1 © © 1 © 2. 

Tensoring further with irreducible representations of the orbital angular momentum op- 
erator L we deduce that for / > 2, 

/ © (0 © 1) © 1 = / © (0 © 1 © 1 © 2) 

= /©(/- 1) ©/©(/ + 1) © (/ - 1) ©/©(/ + 1) 
© (/ - 2) © (/ - 1) © / © (/ + 1) © (/ + 2). (3.11) 

For / = 1 we have 

1 © (0 © 1) © 1 = 1 © (0 © 1 © 1 © 2) 

= 1©0©1©2©0©1©2©1©2©3, (3.12) 

while for / = we have 

0©(0©1)©1 = 0©1©1©2. (3.13) 

We can use these equations to count the number of angular momentum representations 
that can occur for a given value of the total angular momentum j. To do this, we fix a 
value of j and count, with multiplicity, the values of / which occur on the right hand side 



of equations (3.11) - (3.13). For j > 2 we find that / = j contributes four times, I = j — 1 
and I = j + 1 three times each, and I = j — 2 and I = j + 2 once each, giving a total of 
12 modes. For j = 1, the value / = j ' + 2 = 3 does not contribute, and / = j = 2 = — 1 is 
impossible, so only 10 modes occur. Finally, for j = 0, the possibility / = occurs once, 
the possibility / = 1 twice and / = 2 once, giving a total of four modes. For a given value 
of j, each of the modes has the usual, additional degeneracy of 2j + 1. However, because 
of overall invariance of the situation under generalised rotations, these 2j + 1 states are 
physically equivalent (and obey the same differential equation). 
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3.3 The J = sector 

The four modes with j = can be constructed very simply by combining / = 0,1,2 
functions with the quaternionic and isospin degrees of freedom to obtain overall scalars 
under the action of J. Since all states must have isospin 1, the I = (constant) function 
can only be combined with e-t to obtain an overall scalar. Using the cartesian coordinates 
of the unit vector x on the spatial two-sphere as the three I = 1 states, we can obtain one 
overall scalar from the scalar (s = 0) field as x-t. Another overall scalar can be constructed 
from the vector (s = 1) field as x- (e x t). Finally, the five independent functions spanning 
the / = 2 multiplet are the components of the tensor XiXj — ^6ij. We can combine these 
with the isovector and the spin 1 part of the quaternion as (x-t)(x-e) — |e-t to obtain an 
overall scalar. We thus have four basis states of the j = sector, and can use any linear 



combination (with coefficients being functions of r) to study the j = sector of (3.8). It 
turns out that, for our purposes, the following basis states are most convenient: 

t;i = (x-t)(x-e) -e-t (3.14a) 

V2 = (x-t)(x-e) (3.14b) 

t;3 = x-(ext) (3.14c) 

V4 = x-t. (3.14d) 



We note that vs and are perturbations of the hedgehog fields as defined in ( |2.6 ) and 
are the assumed shape of the perturbations used in [15]. 

It is worth interpreting our chosen basis physically before proceeding with the math- 
ematical analysis. The generator x-t in the isospin Lie algebra is in the direction of 
the asymptotic Higgs field and thus the generator of the unbroken f/(l) subgroup of the 
original isospin symmetry. Even though the YMH model studied here is not a realistic 
physical model, we adopt a terminology where this U{1) is interpreted as the gauge group 
of electromagnetism. Then we note that V2 is the unbroken part of the SU2 gauge field 
and can therefore be thought of as the photon field. The other spin 1 states Vi and fa are 
eigenstates of (x-t)^ with eigenvalue —1 and therefore have electric charge ±1. Neither 
of them is an eigenstate of x-t so both should be thought of as linear combinations of 
excitations of the charged ly-bosons in this model. Finally, the spin zero state ^4 is pro- 
portional to the asymptotic value of Higgs field and describes a massless and uncharged 
excitations of the Higgs field. We still need to ascertain which linear combination of these 



states satisfies the background gauge condition (2.26). We will do this below by applying 



the operator Jp, as outlined at the end of Sect. 2.4 



Our next task is the application of the operator {!/) IJ))bps (3.3) to a linear combination 



ci = ^hn{r)vn, (3.15) 



n=l 



of the basis functions of the j = sector. In order to organise the calculation we note 
that {Ip^ 1P)bps is a linear combination of the operators A, L - 1, t^, (x - 1)^, (x - e)(x - 1) 
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Table 1: The action of the operators in (3.3) on each of the zero angular momentum 



modes (3.14). 











Da 




— 2f?;i + 2vr,) 


—2(vi + 2vo) 




—2vi 


L t 


Vl + 2f 2 


Vl + 2f 2 


^^3 


2vi 




-2vi 


-2172 


-2v^i 


-2vi 


(x-t)2 







-V3 





(x-e)(x-t) 


Vl 





V3 





et 


Vl - 2'U2 


-Vl 


V3 - 2V4 


-V3 



and e ■ t, multiplied by simple functions of W and H. We list, in Table 3.3 the action of 



the relevant operators on the modes wi, . . . ,^4 (3.14). 



Using (3.3), we find that (3.8) becomes 



— r{rhi)"vi — r{rh2)"v2 — rirh^y'v^, — r{rb4)"vij + {^biVi + 262^1 + 46if2 

+ Ab2V2 + 2b3V3 + 2641^4) + ( - 2(1 - W)biVi - 2(1 - W)b2Vi 

- 4(1 - W)biV2 - 4(1 - W)b2V2 - 2(1 - W)b3V3 - 4(1 - Vr)M4) 

+ (2(1 - W)\vi + 2(1 - W)%V2 + 2(1 - W)%V3 + 2(1 - W)%4Va^ 

+ ((-(1 - + H')bivi + (-(1 - + H')b3V3) 

+ {2{rH' - (1 + W)H)biVi + 2{rH' - (1 + W)H)b3V^ + {2WHbiVi 

- 2WHb2Vi - 4WHbiV2 + 2WHb3V3 - 2WHbiV3 - mHb3Vi 

Pu'^{biVi + b2V2 + b3V3 + b^Vi). 



(3.16) 



Each group of terms which is bracketed in (3.16 ) comes from one part of (3.8 ). We compare 



coefficients of fi,f2,f3 and V4 in (3.16) to obtain equations for the radial functions 6n(r). 
We find that 61 (r) and b2{r) are coupled and that 63 (r) and 64 (r) are coupled. As noted 
previously, bsi^r) and 64(r) are related to the functions w{r) and h{r) in [T5j, as they are 



also perturbations of the basic hedgehog fields (2.10 ). The coupling is not quite symmetric 
but this can be remedied by defining 



w 



rb3, h 



A similar redefinition 



rbi, 



a 



rb^ 
T2 

rb2 

T2' 



(3.17) 



(3.18) 



for the other coupled system aids comparison with the previous work, and again results 
in the coupling being symmetric. 
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With these abbreviations, we conclude that the insertion of 

1 , 



V4 



(3.19) 



into (3.8) gives a system of second order differential equations which decouples into two 
systems: 



^ 3W^ + 



dr 



V — 



-a 



2W^ + 2\ 2^W{H-l) 



a 



and 



^2 3iy2 ^ JJ2 



dr 



dr'' 



+ 



2W''' 



2^WH , 

w n 



, 2^/2WH 

h w 



(3.20a) 
(3.20b) 

(3.21a) 
(3.21b) 



As noted, w{r) and h{r) are perturbations of the hedgehog fields W{r) and H{r), and 
thus we could have obtained (3.21) by linearising (2.16). This is what is done in [15]. In 
the following sections we will investigate bound and scattering states in the systems (3.20) 
and (3.21). Recalling the interpretation of the basis (3.14) we note that the system (3.20) 
describes a photon mode interacting with a IV-boson mode, and that (3.21) describes a 
massless Higgs perturbation interacting with a W^-boson mode. Correspondingly, we will 
often refer to the former as the photon system and the latter as the Higgs system. 



So far, we have focussed on the equation (2.42) since it is directly related to the linearised 
YMH equations. The alternative but equivalent equation (2.45), however, also plays an 
important role in our analysis. As explained at the end of Sect. 2.4, the background gauge 
condition 2.26 is conveniently implemented in this formulation. It also turns out that the 
equation(2.45) is also sometimes easier to analyse than (2.42). 



Recalling the definition (3.6), we thus consider the eigenvalue problem 



(3.22) 



Inserting, in analogy to (3.19), the expression 

1 



P = -{xvi + V2iJV2 + ^V3 + V2C 



(3.23) 



for four radial functions ^, r/, x, we obtain again two systems of second order differential 
equations: 



d^ W^ + H^ + l 
dr"^ ^2 

rf2 2W^ 
dr"^ 



(3.24a) 
(3.24b) 
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and 



W'' + H'' + l\ 2^fm , 2 
H 1 X H -2 — ^ = ^ X, 



J2 



-r^ + 5 ^ + = ^ ^■ 



(3.25a) 
(3.25b) 



According to (2.46), we can relate a solution (3.19) of (3.8) to a solution (3.23) of (3.22) 
via the Dirac operator JP'^bps (3-1) according to 



(3.26) 



We are now in a position to apply the background gauge condition (2.26) to solutions 



of (2.45). As discussed at the end of Sect. 2.4, we can do this by ensuring that the 



quaternion p has no real quaternionic part. Of all the basis elements (3.14), only is 



real in the quaternionic sense. Hence solutions of the form (3.23) correspond, via (3.26) 



to solutions of the linearised YMH equations if the coefficient function C, in (3.23) is 



identically zero. Note that, according to (3.24), this can be imposed consistently since C, 



satisfies a homogeneous linear equation and does not couple to any other mode. 



Working out the relation (3.26) explicitly is a little tedious, but can be done by careful 



application of some of the results in (3.3) as well as the application of the operator {cidi) 
to fi,...,f4. It can be seen without too much effort that the transformation (3.26) 
permutes the basis elements f i, . . . , t>4. In particular, the coefficient functions ^, C, of ^3, U4 
are mapped to the coefficient functions v,a of Vi,V2- Equally, the coefficient functions 
X, ijj of f 1, f 2 are mapped to the coefficient functions w, h of fs, ^4. Since the background 
gauge condition does not restrict the functions x ^-iid ip, we deduce that any solution of 



(3.21) satisfies the background gauge condition. 



However, in order to understand the implication of the background gauge for (3.20), we 



need to know the relation between ^, and f,a. To compute the effect of applying (3.1) 



to the basis (3.14) we note that, in terms of the usual spherical coordinates {r,9,(p) for 



^sin 6 cos if^ 
sin 9 sin Lp 
cos 9 



9 = de± 



^cos 9 cos (f^ 
cos 9 sin 
— sin 9 




and 



X X L = 9dt) 



sm 9 



Therefore 



Cidi = it-edr + -9-edg H ^^—-ip-ed^ = it-edr + -e-x x L. 

r r sin 9 r 



Then we compute, for example. 
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Evaluating the other terms (3.26) we arrive at 



uja = + ^, 

ar r r 

UJV = 4 C- 

ar r r 



(3.27a) 
(3.27b) 



Remarkably, the transformation (3.27) relates the coupled system (3.20) to the decoupled 



system (3.24). Thus, the most efficient way of studying the the coupled system (3.20) is 



to look at (3.24) instead, and then apply (3.27). In addition, if one is only interested in 
bosonic modes satisfying the background gauge condition (as we are in this paper), one 
should set ( to zero. 



4 Spectral properties of the j = sector 
4.1 General considerations 

Before we turn to the numerical investigation of bound states and scattering in the j = 



sector, we note some general features of the systems (3.20) and (3.21) on the one hand 



and (3.24) and (3.25) on the other. For this purpose we note the asymptotic behaviour 



of the coefficient functions W and H (2.10): 

W{r) ^ 



r 
~6 



H r) 



r 

T 



for small r, 



(4.1) 



and 



W{r) 



H{r) 1 — r + 0(e for larger. 



(4.2) 



It follows that the systems (3.20) and (3.21) both decouple at large values of r. We 



can think of them as two channels which are coupled in the region of the background 
monopole but whose coupling falls off exponentially fast as we go away from the core 
of the monopole. It is instructive to set the coupling terms to zero and consider the 



resulting single channels. The equations (3.20a) for v and (3.21a) for w are the same after 
decoupling, and take the form 



dr"^ 



(4.3) 



The potential appearing in this equation has the asymptotic form 



+ 



l-- + 0{e 



as r — )■ oo, 



(4.4) 



and, in particular, tends to the positive constant 1 for large r. Thus, thinking of (4.3) as 



the spatial part of a radial wave equation, we see that its wave solutions correspond to 
massive particles of mass 1, which is in agreement with our interpretation of v and w as 
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excitations of the VT-bosons in YMH theory. By contrast, the equation for the function 
a after decouphng contains the potential 



+ 22 

~ ^ + Oie for larger, (4.5) 

while the equation for the function h after decoupling contains the potential 

2W^2 

0(e ^) for larger. (4-6) 



Therefore, the corresponding radial waves are massless excitations, in agreement with 
their interpretation as, respectively, photon and Higgs excitations. 

The potential + - l)/r2 is plotted in Fig. |l| The plot and the appearance of 

the attractive Coulomb tail in the asymptotic form ( |4.4[ ) suggest that the Sturm-Liouville 
problem (4.3) should have an infinity of bound states for o;^ < 1 accumulating at o;^ = 1, 
with a continuous spectrum for u"^ > 1. 




The coupling of massless channels to a massive channel with infinitely many bound states 
constitutes the generic situation for the occurrence of Feshbach resonances, or quasi- 
normal modes. We give a brief summary and references in Appendix |Xj On the grounds of 
the general theory, we might expect both (3.20) and (3.21 ) to exhibit resonance scattering 
for < cj^ < 1. However, we also note an important difference between the systems (3.20) 
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and (3.21). It follows from (4.1) that the coupling terms in (3.20) are singular at r = 



whereas they are smooth for (3.21). 



Turning now to the systems (3.24) and (3.25) we observe they, too, consist of a massive 



and a massless channel each. The potential in the massive channels is 



2 2 

- + 4 + 0(e- 



as r — 7- oo, 



(4.7) 



and, because of the attractive Coulomb term, should again support infinitely many bound 



states. However, while the system (3.25) has the same generic form as the photon and 



Higgs systems (3.20) and (3.21) discussed above, the equations in (3.24) are decoupled. 



As a result, we expect this system to have true bound states for the excitation denoted 
^ in the range < < 1. There are scattering states in the excitation (, but these 



do not satisfy the background gauge condition, as already discussed after (3.26). Hence 



on the basis of the decoupled form (3.24), we expect the system (3.20) to have infinitely 
many bound states and no (physical) scattering states in the range < < 1. This is 



in marked contrast with the Feshbach resonances of the Higgs system (3.21) 



Before we present our numerical results we briefly comment on or numerical methods. 
All of our scattering calculations require a shooting-to-a-fitting-point method for finding 
solutions which decay at large enough values in one of the channels. Numerical instability 
would always create difficulties where the exponentially increasing term creeps into the 
solution. Thus, we choose an appropriately large value of r to integrate to and use a 
shooting method to find the decaying solution. 

4.2 Bound states 

It is well known that the zero-modes of the 't Hooft-Polyakov monopole give rise to zero- 
energy bound states of the linearised YMH equations. The zero-modes of the 't Hooft- 
Polyakov monopole are obtained from infinitesimal translations in and a special 'large' 
gauge transformation generated by the Higgs field itself. This transformation does not 
vanish at infinity and is therefore considered as a physical symmetry transformation rather 
than a gauge transformation. The infinitesimal form of the 'large gauge transformation' 
generated by the Higgs field is 



0, 6Ai = D4 = -B,. 



In our quaternionic formulation the corresponding zero-mode is simply 



while the other three zero-modes related to translations are 



(4.9) 



ik = ek{eiBi 



(4.10) 
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It is easy to check that (4.9 ) (and hence (4.10 )) satisfy the hnearised Bogomol'nyi equation 
(2.35), and hence (3.8) for u = 0. Of the four zero-modes just found, only has j = 0. 



It has the explicit form 



^ ^ , XiXa XiXa , ^^^-.H WH 
Bi = -A0 = -^H'ta + ^1 + W)^ta ^U, 



(4.11) 



and is thus a linear combination of the basis states vi and V2 as defined in (3.14). In 



particular, it is therefore a solution of the photon system (3.21). 

Bais and Troost first investigated bound states in single channels arising in the linearised 
BPS system in |8j. In particular, they computed bound state energies in the system 



(3.24a). For completeness of our discussion we have repeated the numerical analysis here. 
Using the NAG shooting method D02KEF for Sturm-Liouville type problems we find the 
eigenvalues and hence ujn > of the first few bound states of (3.24a), to 3 significant 



figures. As explained in our qualitative discussion in the previous section, we expect 
there to be infinitely many Coulomb bound states in this channel. We can estimate 



their energies by neglecting the exponentially small terms in the potential (4.7) and using 



the standard formula for Coulomb bound state energies. In Table 4.2 we list both the 



numerically computed bound state energies of (3.24a) and the Coulomb approximations 



1 



(4.12) 



(n + l)2' 

where we used the standard expression for Coulomb bound state energies, recalling that 



/ = 1 in this case by virtue of the 1/r^-term in (4.7). Even for very small ra, the Coulomb 
approximation is surprisingly good. 



n 




1 


2 


3 


4 
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0. 


877 


0.946 


0.970 


0.9^ 


^0 


0.997 




0. 


866 


0.943 


0.968 


0.9^ 


^0 


0.996 



Table 2: Values for Un for the eigenvalue problem (3.24a), and Vtn given in (4.12) 



We have also computed the wavefunctions 



n 



1,2,... for each of the bound states. 



They have the standard form of Coulomb bound states, but it is interesting to note that 



they correspond, via (3.27), to bound states in the coupled system (3.20). The functions v 



and a for each of the bound states describe, via (3.19), the profile of the excited monopole. 
Since h and w both vanish for this excitation the bound state only involves the photon 
excitation a and the component v of the field describing the W-boson. 
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From the point of view of the system (3.20), it is surprising that there are true bound 
states. As explained in Sect. 4J^, one would generically expect two-channel problems like 
(3.20) to have Feshbach resonances, but no bound states. 



4.3 Scattering 

In Appendix |A] we discuss what kind of coupled systems are likely to possess Feshbach 



resonances. Both (3.20) and (3.21) have the required features in the parameter range 
< cj^ < 1: they consist of two channels which are weakly coupled at large distances, and 
are such that, after removing the coupling term, there are bound states in one channel 



and only scattering states in the other. In fact, the system (3.21) shows exactly the 



expected Feshbach resonance behaviour. This was first noticed by Forgacs and Volkov in 



[I5] . and we will revisit it below. However, we will also see that the system (3.20), whose 
bound states we analysed in the previous section, has no scattering states which satisfy 
the background gauge condition. 

We begin with a brief review of the results of [15] . This will establish our conventions and 



terminology and serve as a preparation of our generalisation in Sect. [5j The system (3.21 ) 
has a regular singular point at the origin, which means that the numerical integration 
must start a little away from r = 0. We make a series expansion near r = and find the 
following leading terms 

w ~ Ar'^, h ^ Br^ for small r, (4-13) 

with real constants A, B. Because of the linearity of the problem we can scale the solution 
to fix one of those constants. The remaining free constant plays the role of the shooting 
parameter. 

Once the correct initial conditions are determined to ensure a decaying solution w (for 



< 1), we can consider the scattering problems in the massless channel (3.21b). It has 
the general form 

£ l{l + 1 
dr"^ r 



2 + + V^(r) - ) u(r) = 0, (4.14) 



where ^(r) — )■ as r — t- oo exponentially fast. In the equation for h, the total potential 
is (|46| so / = 0. 



For large values of r, where V can be neglected, the solution must be a combination of 
spherical Bessel functions 

u{r) = r {Aiji{ur) + Biyi{ur)) , for r large. (4-15) 

Since, asymptotically, ji{cor) ^ sin(u;r — ^7r)/{ujr) and yi{cor) ^ — cos{ujr — ^7r)/(ujr) we 
have the asymptotic form 

^ ^ ^Af + Bf . ^ , ^ ^ 

u{r) ^ — — — sm{ujr + O; — -tt) 

2 
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of the radial wavefunction, with the phase shift defined via 



tan((5i 



By evaluating both sides of ( 4.15[ ) for two large values of r (or by evaluating both sides 
and their derivatives for large r) we extract the coefficients Ai and i?;, for each value of 
uj , and hence the phase shift 5i dX uj. The latter determines the partial scattering cross 
section according to the standard expression 



47r(2/ + 1) 



sin 5i{uj) 



and the total cross section according to 



An 



00^ 



^(2/ + l)sin2 5i{uj). 



1=0 



Near a resonance, 6i increases rapidly by it. The function sin^ 6i{u) takes values between 
and 1, is maximal at the resonance and is thus an expedient quantity to plot when 
looking for resonances. 

Applying this procedure to the / = contribution to the scattering cross section for the 



massless channel h from the Higgs system (3.21) we confirm the result shown in Fig.l 



of [T3], suggesting infinitely many resonances as the energy approaches the critical value 



1. A graph of sin^ 5(ci;) against u for the massless channel h is shown in Fig. [2] 



We have also computed the bound state energies of the decoupled massive channel (4.3) 



which occurs in both (3.20) and (3.21), using the numerical method summarised in the 



previous section. They were also computed in [15]. Our results are listed in Table 4.3 



given to 4 significant figures. They are in good agreement with the results in Table 11 of 
[T^ . Comparing with Fig.|2| one sees that the energies of the bound states are close to the 
energy values where the resonances occur. The bound states of the decoupled problem 
have turned into resonances in the coupled problem. This is typical Feshbach behaviour. 



n 1 


2 3 4 5 10 
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Un 0.7984 


0.9263 0.9618 0.9766 0.9842 0.9956 


0.9980 


Table 3: 


Values of Un for the eigenvalue problem ( 


4.3 


)• 



Next we turn to scattering states in the system (3.20 ), exploiting the equivalence with the 
simpler system (3.24). In the energy range that we are interested in, with u"^ < 1, only the 
equation (3.24b) has scattering solutions. However, from our discussion after (3.26) we 
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0.8 0.85 



Figure 2: For the system (3.21), the partial scattering cross section sin^ 6o{uj) is plotted 
against u. 



know that a non-vanishing ( violates the background gauge condition. We can therefore 



conclude that the system (3.20) has no scattering states satisfying the background gauge 
condition in the range < 1. This is surprising in view of the superficial similarity to 



the system (3.21), which, as we saw above, shows interesting resonance scattering. 



We end this section with some observations about scattering solutions of (3.24b). Even 



though they do not satisfy the background gauge condition and therefore do not corre- 
spond to solutions of the linearised YMH equations we expect them to play a role in 
the supersymmetric version of the theory as fermonic scattering states. The scattering 



problem associated to (3.24b) is also of independent interest since it can be solved exactly. 



Inserting the expression (2.10) for the profile function W, the equation (3.24b) takes the 
following simple form: 

2 



+ 



sinh^( 



-c 



(4.16) 



This equation was studied in the context of fermion scattering off monopoles in (the 



authors there did not include the Higgs field, but this does enter (4.16) anyway). A 
scattering solution regular at the origin is 



C(r) = {iu + coth(r))e-*'^'^ + {zoo - coth(r))e*'^", 
from which we can read off the expression 



i + uj 



I — u 



(4.17) 



(4.18) 



for the phase shift. Hence siv? 5{uj) = ^^/(l + u"^), so that the partial scattering cross 
section is simply 



An 
l + w2 
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5 Perturbations of the non-BPS monopole 



We now extend our study of the linearisation around SU (2) monopoles in YMH theory by 



moving away from the BPS hmit and allowing for a non-zero value of A in (2.3). Outside 



the BPS limit we cannot use the quaternionic version of the linearised YMH equations, 



(2.36), because it relied on the quaternionic version of the Bogomol'nyi equation, (2.35) 



which no longer holds. We therefore will not study the most general perturbations in this 



section but restrict attention to the hedgehog ansatz (2.6) and study perturbations within 



this ansatz. With A 7^ we do not have the analytic solutions (2.10) for W{r) and H(r 



of the equations (2.7) for the profile functions in the Hedgehog, so we must solve these 
numerically. We use the NAG routine D02GAF, which solves non-linear boundary value 
problems using a finite difference technique with deferred correction. Plots for a range of 
different values of A can be found in 

As we switch on A, the asymptotic behaviour of the monopole changes, and this will be 



crucial in studying perturbations. The asymptotic function for W{r) is still as in (4.1), 
but H{r) no longer has a constant term as r — 00. The new leading terms are 

W{r)^^e~'', H{r) ^ -r + 0{e~'^) as r 00. (5.1) 

Even for small values of A, integration beyond around r = 10 becomes unstable, so for use 
in the linearised equations we spline with the asymptotic approximations for W{r) and 
H{r). In studying the linearised problem we fix a value of A and use the profile functions 
for that case. Thus, from now on we set A = 0.1. 



The linearisation of (2.7) was already given in [15]. Denoting the hedgehog profile func- 



tions by and W^, inserting 

W{r, t) = W'{r) + e^"*w(r), H{r, t) = H' + y/2e'''%{r) (5.2) 
and keeping only linear terms in w and one obtains, after re-naming — H , — )■ W , 



H w — ly 1 — —h = bj w, 



d^h 2W^ , , 



dr 



(5.3a) 
(5.3b) 



We note that (5.3a) is the same equation as (3.21a), but the differing asymptotic behaviour 
(5.1 ) of H{r) means that the potential appearing in this equation now behaves as 



~ 1 + Die ) tor large r. 



(5.4) 



This potential is plotted in Fig. [3] 
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Figure 3: The potential ^ plotted against r, for A = 0.1. 



The system (5.3) has the same generic form as the systems (3.21) and (3.20) studied in 



the previous section. We follow the same strategy in analysing it as we did in Sect. 4.3 



Thus we first consider the bound state problem, where the right hand side of (5.3a) is set 



to zero, so that it is artificially decoupled from the system. The eigenvalue problem we 
need to solve is 

(5.5) 



dr' 



+ 



w 



As we saw in (5.4), the potential appearing here tends to the limiting value 1 faster than 



in the Coulomb case but still according to an inverse square law, which is the threshold 
case for supporting infinitely many bound states, as discussed in in [21]. As shown there, 
the radial Schrodinger problem 

d^w 



dr' 



+ V{r)w = Ew, 



(5.6) 



where V{r) ~ ^ for large r, has infinitely many bound states with i?<0if/3<— | 
(for > /3 > — I only a finite number of bound states are supported). Thus we expect 



infinitely many bound states with w < 1 in (5.5), with all but the lowest values of oj 



lying very close to 1. This poses problems for the NAG routine D02KEF, which requires 
a good initial estimate of the eigenvalue. Fortunately, the inverse square potential is 
a well-studied problem, and the WKB method provides an approximation for the ratio 



of successive eigenvalues [22]. Writing u"^ for the n-th eigenvalue in (5.5) we have the 
following formula, valid for large n: 

1 



OJ. 



n+l 



1 



47r 



(5.7) 



The first eigenvalue and eigenfunction is easily computed using the NAG routine D02KEF; 
the eigenfunction is displayed in Fig. |4} Initial guesses for subsequent eigenvalues can then 
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be computed using (5.7). This approximation expected to be good only for large ra, but 
worth using as an initial estimate in D02KEF for n = 2, 3, 4. The computed eigenvalues, 
predicted values, and estimated error in computed value are displayed in Table |5| Rather 
surprisingly, the WKB guess even for = 3 is already accurate to 10 decimal places. The 
associated eigenfunctions are plotted in Fig. |5| They only show the exponential decay 
characteristic of bound state wavefunctions for rather large values of r, and we therefore 
plot them twice, using different scales. 




Figure 4: The lowest bound state of (5.5) for A = 0.1 



n 


1 


2 


3 




4 




0.95343 


0.99998 


0.9999999< 


35 


0.999999999990 


WKB prediction 


n/a 


0.99997 


0.9999999< 


B5 


0.999999999989 


duJn 


0.38 X 10-6 


0.68 X 10-6 


0.63 X 10- 


10 


0.42 X 10-13 



Table 4: Eigenvalues Un of (5.5) with A = 0.1, the predicted values for n > 1 from 5.7 



and the estimated accuracy duJn of the computed value a;„. 



As discussed in appendix A the presence of bound states in (5.5) suggests the presence of 



(Feshbach) resonances in the radiative channel of the two channel problem (5.3). We now 



look for such resonances by studying the scattering associated with (5.3), in the energy 
range < 1. We follow the procedure summarised and used in Sect. 4^ Thus we 
integrate (5.3) for A = 0.1 and a range of u"^ < 1, using the D02L suite of NAG routines 



for the integration. These solve initial value problems for ODEs using the Runge-Kutta- 
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Figure 5: The 2nd (top) and third (bottom) bound state of (5.5), for A = 0.1. The scale 



on the right is chosen to show the first oscillation near the origin 



Nystrom method. We then tune the initial conditions with a shooting method to ensure 
that the function w{r) decays at large r. The system has a regular singular point at the 
origin, which means that the numerical integration must start a little away from r = 0. 
We make a series expansion of (5.3) about r = to find the correct initial data. We find 
that w{r) ^ Ar"^ and h{r) ^ Br^ for real constants A, B, as in the A = case. 

We determine the ratio A/B which allows w to decay at infinity and integrate to compare 



the solution for h at large distances to the spherical Bessel function Jq, as in (4.15) and 



to extract the phase shift 6{uj). Our results are shown in Fig. [6j We see that as a; — )■ 1, 
the phase shift 6 increases sharply in value. An increase in 5 by tt suggests a resonance. 
However, successive and closely spaced increases cannot easily be separated. 

In Fig. [t] we plot sin^ 5 as a function of co. We can clearly see one fairly wide resonance 
centred around u = 0.985. Again, the accumulation of the bound state energies at 1, as 
shown in Table |5| makes it difficult to discern the interesting behaviour as w — 1. We 
stretch out the region near threshold to see if there are, as we expect, more resonances 
as w — > 1. The result is plotted on the right in Fig. [7) We can see that there is at least 
one more resonant peak. We solved (5.3) for energies up to a; = 0.999998, which takes us 



past only the first two bound state energies in Table [5} thus the occurrence of two distinct 
resonance peaks in the scattering cross section is exactly what is expected. At higher 
energies the wavefunctions are extremely long range, of the order of r = 10^. Therefore, 
integrating for higher energies becomes increasingly time consuming. We expect, however. 



that each of infinitely many bound states of (5.5) will turn into resonances, and that 



infinitely many more resonant peaks in the scattering cross section arise as w — )■ 1. 
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6 Discussion and conclusion 



The linearised YMH equations in the background of a 't Hooft-Polyakov monopole have 
very interesting spectral properties. We have seen that there is a zero-energy bound 
state and an infinity of Coulomb bound states in the same energy region as an infinity 
of Feshbach resonances. For vanishing total angular momentum, there are two coupled 
systems which, despite their superficial similarity, have very different spectral properties: 
in one (called the Higgs system here) there are infinitely many Feshbach resonances, while 
in the other (called the photon system here) there are infinitely many true bound states, 
but no scattering states satisfying the background gauge condition and, in particular, no 
resonances. The occurrence of bound states can be understood in terms of a decoupling 
transformation. We also showed that the Feshbach resonances in the Higgs system persist 
when the Higgs self-coupling A is switched on. We saw that they become even more densely 
spaced, essentially because an attractive 1/r potential in the BPS limit is replaced by an 
attractive potential when A 7^ 0. 

In this paper we restricted attention to the frequency range u"^ < 1 and to spherically 
symmetric perturbations (in the generalised sense). In this setting, the scattering is 
effectively single channel scattering, fully described by a single phase shift. When one 
goes beyond the threshold u"^ = 1, the scattering will be genuine two-channel scattering, 
whose study is more involved, both numerically and in terms of the interpretation of 
the results. However, very similar multi-channel scattering problems are much studied 
in atomic and nuclear physics, and were considered in [23] in the context of monopole 
scattering, so there is no problem in principle of carrying out a similar study here. 

It would also be interesting to explore systems arising for larger eigenvalues of the gen- 
eralised total angular momentum operator J. The combination of our quaternionic for- 
malism with the techniques developed in [5l El El] should provide an efficient method 
for finding the corresponding systems of coupled differential equations. As explained in 
Sect. |3.2[ we expect the system for j = 1 to consist of ten equations, and for j > 2 of 
twelve. However, a parity argument will split each of these system into two sub-systems 
(as it did in our j = case), so that the largest system one needs to consider has six 
channels. 

To end, we point out the striking similarity between the spectral properties of the lin- 
earised YMH system in the background of the 't Hooft-Polyakov monopole and those 
of the Laplace operator on the moduli space of charge two SU{2) monopoles [23l 120] . 
The latter also include zero-energy bound states. Coulomb bound states embedded in the 
continuum and Feshbach resonance scattering. These similarities are likely to have an 
interpretation in terms of electric-magnetic duality. Spelling this out is left as the topic 
for a future investigation. 
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A Feshbach resonances 



A Feshbach resonance ^3] is a resonance in a system consisting of several channels in 
which a bound state occurs if the coupling(s) between the channels vanishes. Feshbach 
resonances are much studied in the context atomic physics [26] but also occur in other 
contexts, see [27] for a pedagogical and recent account. 

Consider a simple system consisting of two channels with coupling terms including, for 
convenience, a parameter < g < 1, 

1 

r^iru) + V(r)u + qC(r)v = Eu, 

r dr'' 

1 d"^ 

—{rv) + V{r)v + qC{r)u = Ev. (A.l) 

r dr'^ 

We can decouple the equations by setting the parameter q to zero. With < 0, suppose 
that the eigenvalue problem 

1 (i^ 

r^(ru) + V(r)u = Eu (A.2) 

has bound states for E < 0, occurring at Eo, Ei, E2, . . . but that those values are part of 
the continuous spectrum for the other decoupled equation 

r dr 



^^{rv) + V{r)v = Ev. (A.3) 



When the coupling parameter q is non-zero, the bound states in (A.2) can leak into 



the radiative channel (A.3) and decay. When that happens, Feshbach resonances in 



(A.l) generically occur at values close to £'0, Ei, E2, . . .. Examples of such resonances are 



studied in [15] and [23], both in the context of magnetic monopoles. A single channel 



eigenvalue problem such as (A.2) is generally much easier to solve than a two channel 
coupled problem such as (A.l). Hence, it is worth looking at the bound state problem, 
as it gives an idea of where the resonances will occur. In systems where the coupling 
terms are weaker at large distances than the potential, the decoupled problem is a good 
approximation. Then these preliminary calculations mean that the region of searching for 
resonances is narrowed down, so that computational time is reduced. 
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